function getEdgeGrid( self )


%Nodal Grids
nodeX = zeros(self.ncx+1,1);
for ii = 1:self.ncx
    nodeX(ii+1) = nodeX(ii) + self.hx(ii);
end

nodeY = zeros(self.ncy+1,1);
for ii = 1:self.ncy
    nodeY(ii+1) = nodeY(ii) + self.hy(ii);
end

nodeZ = zeros(self.ncz+1,1);
for ii = 1:self.ncz
    nodeZ(ii+1) = nodeZ(ii) + self.hz(ii);
end

%Cell centre grids
ccX = nodeX(1:end-1) + self.hx/2;
ccY = nodeY(1:end-1) + self.hy/2;
ccZ = nodeZ(1:end-1) + self.hz/2;

[Xex, Xey, Xez] = ndgrid(ccX,nodeY,nodeZ);
[Yex, Yey, Yez] = ndgrid(nodeX,ccY,nodeZ);
[Zex, Zey, Zez] = ndgrid(nodeX,nodeY,ccZ);

self.Xe = [Xex(:), Xey(:), Xez(:)];
self.Ye = [Yex(:), Yey(:), Yez(:)];
self.Ze = [Zex(:), Zey(:), Zez(:)];

end

